The Equilibrium Intrinsic Crystal-Liquid Interface of Colloids 
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We use confocal microscopy to study an equilibrated crystal-liquid interface in a colloidal sus- 
pension. Capillary waves roughen the surface, but locally the intrinsic interface is sharply defined. 
We use local measurements of the structure and dynamics to characterize the intrinsic interface, 
and different measurements find slightly different widths of this interface. In terms of the particle 
diameter d, this width is either 1.5d (based on structural information) or 2.4d (based on dynamics), 
both not much larger than the particle size. This work is the first direct experimental visualization 
of an equilibrated crystal-liquid interface. 

PACS numbers: 



INTRODUCTION 

The interface between crystal and liquid phases of a 
material governs phenomena such as wetting, lubrica- 
tion, and crystal nucleation [l], Q • Interfaces are poorly 
defined at the atomic level: capillary waves cause fluc- 
tuations in the interface position 0, 0, S [B| , and locally 
the structure varies in a smooth way from ordered to dis- 
ordered Q. The literature makes a distinction between 
the intrinsic interface (presumed to be sharp) [3j , and the 
observed surface blurred by capillary waves [1,01. The 
standard equilibrium interface profile is well-defined and 
contains fluctuations at all length and time scales, due to 
these capillary waves. Because of the rapid time scales of 
capillary wave fluctuations, along with the small length 
scales at the interface, it is difficult to study these in- 
terfaces directly [lj. Thus computer simulations provide 
useful information about model crystal/liquid interfaces, 
such as hard sphere systems [3, la, IS) an d Lennard- Jones 
systems [13, Ell- 

Recently, crystal/liquid interfaces were directly studied 
in colloidal suspensions using confocal microscopy [rl[l2|. 
Colloids are systems of solid particles in a liquid, and are 
a good model system for phase transitions [H, 0, [IH, HH . 
Microscopy allows direct observation of structure and dy- 
namics of the colloidal particles [l^|. However, the previ- 
ous experiments focused on nonequihbrium cases where 
samples were crystallizing, and did not provide data on 
equilibrium interfaces such as those studied by simula- 
tion @, II, S 113, EH] . Furthermore, these experiments did 
not examine the intrinsic interface, perhaps because they 
were nonequihbrium studies and thus crystalline particles 
were present in the "liquid" side and vice-versa, which 
confuse the structure near the interface. 

In this work, we present confocal microscope observa- 
tions of an equilibrated colloidal crystal/liquid interface. 
By following the positions of several thousand colloidal 
particles on both sides of the interface, we directly vi- 
sualize the interface. An example of our data is shown 
in Fig. HJa), where blue particles are crystalline and yel- 



low/red particles are liquid-like. This interface has a low 
surface tension and we see capillary waves. We are able 
to remove the influence of these capillary waves from the 
data and measure the intrinsic surface profile. In par- 
ticular, we find that capillary waves cause an apparent 
broadening of the surface, but the structure of the in- 
trinsic surface is characterized by a width of only 1.5d 
(in terms of the particle diameter d). This is the first 
direct experimental visualization of an equilibrated in- 
terface. A precise definition of the intrinsic interface has 
not been made before, and we show that several plausible 
definitions give slightly different results. 

Our colloidal sample is described in the Materials and 
Methods section. Here, we briefly note that in our sol- 
vent the colloidal particles have a slight charge. Experi- 
mentally we observe that the freezing transition volume 
fraction is at </>f ree ze = 0.43 and the melting transition is 
at t^meit = 0.49, to be compared with the hard sphere val- 
ues of ^freeze = 0.494 and </> mo it = 0.545. Our measured 
values are similar to that seen for other experiments with 
similar colloidal samples [l2j |. 



RESULTS 

At each time step, we determine the crystalline re- 
don by using the method of bond-order parameters 
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To do this, each particle i at each time 
is characterized by a normalized order parameter q ; (i) 
with (21 + 1) complex components 



qim{i) = 



NB{i 



B(i) 



(i) 



where M is a normalization factor such that 
J2 m 9im(i)qt m (i) = 1, B(i) is the number of neigh- 
bors of particle i, fij is the unit vector pointing from 
particle i to its j'th neighbor, and Yi m is a spherical 
harmonic function. Following prior work, we use I = 6 
[13, EH- The neighbors of a particle are defined as 
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FIG. 1: (Color) (a) Picture of the liquid/crystal interface, 
showing a slice that is two crystalline layers thick (4.8 /j,m). 
Particles are colored according to the number of ordered 
neighbors they have, N . Blue particles have N > 8 and 
darker blue indicates higher N . Yellow, orange, and red par- 
ticles have N < 8, with red particles having N — and the 
brightest yellow particles having N — 7. The image has been 
slightly rotated (by 3° around an axis parallel to gravity) to 
view the crystalline particles in-plane, (b) Here the parti- 
cles are colored according to how much time they spend as 
"crystalline" particles (with N > 8). White particles spend 
all of their time as crystalline or liquid-like, and the colored 
particles fluctuate over the ~ 1 hr experiment. Blue particles 
spend nearly all their time as crystalline, light green corre- 
sponds to those which spend half their time in each state, 
and red particles spend nearly all their time as liquid-like. 
While the coloring is based on the time average of the data, 
the positions are considered at the same time as shown in 
panel (a). 



those with centers separated by less than lAld, which 
corresponds to the first minimum of the pair correlation 
function g(r) for the liquid region. Two neighboring 
particles are termed "ordered neighbors" if the complex 
inner product J2 m Qim(i)Q* m (J) exceeds a threshold value 
of 0.5. For each particle, the number of ordered neigh- 



bors N Q is determined. Following the usual convention, 
particles with N a > 8 are classed as crystalline par ticles, 
and the other particles are liquid- like particles [lq . The 
advantages of using bond-order parameters is that they 
are local measures of order, are somewhat insensitive to 
variations in the number of neighbors each particle has, 
and also do not depend on the specific type of crystal 
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Figure QJa) shows a snapshot of the sample, colored 
according to the bond-order parameter method. The 
crystalline side of the interface (blue particles) is com- 
posed of hexagonal layers in random stacking (a mixture 
of face-centered-cubic stacking [abcabc...] and hexagonal- 
close-packed stacking [ababab...]). This is similar to the 
stacking seen in growing colloidal crystallites In 
particular, the data shown in this manuscript are for a 
sample where 9 hexagonal layers are imaged, with stack- 
ing abacbebac. Two other regions of this sample were 
imaged and the results presented below do not vary in 
any significant way. 

Typically, between two different phases, there will al- 
ways exist a surface energy defining the energetic cost 
of maintaining the interface. For fluid/fluid interfaces, 
this is also known as the surface tension, although the 
concept is relevant for solid/liquid interfaces as well (for 
example, in nucleation of crystals [l2|.) For a hard sphere 
system, this surface energy is entropic in origin, reflecting 
the difficulty of packing particles at the interface. Right 
at the interface, particles can pack neither optimally for 
the crystal state, nor optimally for the liquid state, thus 
resulting in an entropic penalty that gives rise to a sur- 
face tension. In our experiment, the interface shown in 
Fig. Q] is rough and varies in time due to surface capil- 
lary waves. Because the surface energy is not extremely 
large compared to ksT (the thermal energy, based on the 
temperature T and Boltzmann's constant fog), thermal 
fluctuations allow these capillary waves to be observed in 
our experiment. The temporal variability of the sample 
is indicated in Fig. Q|b) , which colors only the particles 
which we observe to spend some time in both crystalline 
and liquid structures during the duration of our exper- 
iment (3750 s). As can be seen, the exact position of 
the interface thus fluctuates by as much as ~ 4 particle 
diameters (in a "peak-to-peak" sense). 

To determine the surface tension, we follow the pro- 
cedure of Ref. [Ej]. Due to a slight purposeful density 
mismatch (described in the Methods and Materials sec- 
tion), the crystal/liquid interface is nearly parallel to the 
xy plane; we rotate the data by 15° to align the interface 
with the xy plane. Based on the new z axis, at each time 
we locate the maximum z position of crystalline particles 
as a function of x and y (coarse-grained in x and y over 
a distance 1.035c?, the spacing between crystalline layers, 
where d is the diameter of an individual particle). This 
gives us the interface height h(x,y,t). While our data is 
only over a short spatial extent in x and y (9d x 24d), we 
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FIG. 2: The circles show the temporal correlation function of 
interface position, averaged over x and y (Eqn. [2}. The line 
is a fit to Eqn. with interfacial stiffness 7 = 0.93 nN/m 
= 1.2fcsT/d 2 and a capillary time r = 750 s. The small 
k modes give rise to the very steep drop of the correlation 
at short time scales, and the exact shape of the theoretical 
curve for At < 25 s depends on how the numeric integration 
of Eqn. [5] is done. 



have good temporal data (0 < t < 3750 s) and thus we 
calculate the temporal dynamical correlation function: 



g h (At) = ([h'(x,y,t)}[h'(x,y,t + At)}) 



x,y,t 



(2) 



where h'(x,y,t) — h(x,y,t) — (h(x,y,t))t represents the 
fluctuations of the interface about its time-averaged po- 
sition, and the angle brackets {...) x ,y.t indicate an aver- 
age over space and time. We plot this dynamic corre- 
lation function in Fig. [2] (circles). The intercept is at 
gh(At = 0) « l.ld 2 ; thus the root-mean-square width 
of the capillary wave fluctuations is approximately d, as 
compared to the slightly larger amplitude of the peak-to- 
peak fluctuations seen in Fig. [T^b). 

The capillary waves are limited by the interfacial stiff- 
ness 7, rather than the surface tension 7 [lfl]. For a 
crystal / liquid interface, 7 is usually anisotropic and de- 
pends on the crystal orientation. Fluctuations in the sur- 
face thus depend both on the surface tension 7 (related 
to interfacial stretching) as well as second derivatives of 
7 with respect to angles away from the interface normal 
(related to interfacial bending) [H, [3, The orien- 

tational average of 7 is the bulk surface tension 7, but 
measu ring this requires vastly more data than we have 
@, [3, EiJ- To extract the interfacial stiffness 7, we fit 
gh(At) using the results of capillary wave theory [EL [20j] . 
Overdamped capillary waves with wave number k should 
decay as 



exp[-t(jk + gA Pi /k)/(2n)], 



(3) 



with gravitational acceleration g, density difference Api 
(across the interface), and viscosity 77 (equal to the sum 
of the viscosities of the two phases) 20] . By equiparti- 
tion, the amplitude of the Fourier component hk of the 
interface displacement contributes as 



<IM 2 > 



with L the lateral system size, and using the capillary 
length £ = y/j/(gA_pi) 21]. As in Ref. 0], we define 
the nondimensional k = £fc, and then combine these two 
results to calculate the theoretical dynamic correlation 
function as an integral 



9h(At) 



k B T 
2777 



■ rf --e X p[-(fe + ^)AV(2r ) ] 
1 ~\~ k 



7L 2 k 2 + £- 



(4) 



with r = £77/7 as the capillary time, the characteristic 
time scale for decay of interfacial fluctuations [2l]]. I n 
particular, r corresponds to the time scale for the slow- 
est decaying wavelength, which is at the capillary length 
scale £. 

We vary r and 7 to find the best fit to our data, and 
plot this as the solid line in Fig. [2j finding excellent agree- 
ment using 7 = l.2k B T/d 2 . Our value of 7 = 1.20 ±0.05 
(in units of k B T/d 2 ) is a similar order of magnitude to 
previously found values of 7 for hard spheres, which range 
from 0.11 to 0.78 @, % [ij, H Our measured 
7 is larger than these values of 7, perhaps due to the 
above-noted anisotropy of 7, for which we may have a 
stiff direction [§] . 

We note several limitations with our measurement of 
the interfacial stiffness. First, our result is calculated 
for one particular interface orientation, rather than av- 
eraging over many different crystal orientations. Second, 
the theory of Refs. [H, is derived for liquid-gas inter- 
faces, where the surface tension is isotropic. As noted 
above, the local interfacial stiffness is an adequate re- 
placement for the isotropic surface tension Q. Third, 
the theory also assumes that both phases have a well- 
defined viscosity, whereas one of our phases is a colloidal 
crystal. The crystalline phase could be considered as a 
liquid with a very large viscosity, but this would then 
predict a very large capillary time, in contrast with our 
observed r = 750 s. This value is consistent with a vis- 
cosity 77 = 50 mPa-sw 20r/o (compared with the solvent 
viscosity 370). This is approximately half of what could 
be expected for high volume fraction colloidal samples at 
volume fractions <fi — 0.43 — 0.49 (2||, given that 77 here 
represents the sum of the viscosities of the two phases. 
To summarize, our crystal-liquid interface has capillary 
waves which fluctuate in a way consistent with predic- 
tions for liquid-gas capillary waves, but interpretations 
of this observation should be done cautiously. 

Additionally we note that the space-averaged height 
(h(x,y,t)} XiV fluctuates but does not monotonically in- 
crease or decrease, confirming that we observe an equi- 
librium interface rather than a system which is crystalliz- 
ing or melting. However, as can be seen in Fig. [TJb), the 
exact interface position fluctuates within a certain range 
over the duration of our experiment, precluding us from 
an ensemble- averaged sort of measurement. For a longer 
duration experiment, the time- averaged shaded region in 
Fig. []Jb) would presumably meander less and be more 
uniform across the y direction. 
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FIG. 3: Two views of the interfacial profile in the "laboratory" 
reference frame, (a) Interfacial profile, showing the number 
density profile d 3 n(z) as a function of z/d. The smooth curve 
is a hyperbolic tangent fit to the envelope, with 10-90 width 
of 9.6d. The inset shows the crystal orientation used to deter- 
mine this curve, after the sample has been rotated slightly (see 
text for details), (b) Bond-order profile, showing the average 
number of ordered neighbors particles have, as a function of 
z/d. The dotted curve is again a hyperbolic tangent fit, with 
a 10-90 width of 6.08d. By definition, (N ) = 8 at z = 0. For 
both (a) and (b), the data are time averaged over 2000 s (80 
3D images). 



We now turn to the structural details of the crystal- 
liquid interface. As mentioned above, we rotate the data 
around the x axis by 15° clockwise (as seen from the view 
in Fig. [1}, so that the crystal structure is aligned with 
the new z axis; see inset of Fig. EJa). We then plot the 
number density d 3 n(z) as a function of z/d in Fig. E^a), 
and see oscillations similar to what has been seen in ex- 
periments @ and simulations @, [13, EH- The smooth 
line in Fig. Ela) is a hyperbolic tangent fit to the data. 
Following previous work, we define the interfacial width 
as the 10-90 width, where the hyperbolic tangent func- 
tion goes from 10% of its value to 90% of its value. We 
find PF10-9Q — 9.6d, quite similar to what has been seen 
previously @, @| • Some of this width is due to capillary 
waves, and some may be due to particles in the crystalline 
lattice (but near the interface) having larger fluctuations 
around their mean positions |2j ■ We checked for this lat- 
ter influence by low-pass filtering Q the density profile 
of Fig. E£a) and found a similar VFio_g width, although 
it depends on the details of the filter; this filtering also 
shifts the halfway point of the hyperbolic tangent fit to- 
ward the liquid side by Ad. We also calculate the width 
based on the structural ordering, by plotting the average 
number of ordered neighbors {N a ) as a function of z' / d 
in Fig.EJb). Based on this measure, the profile is slightly 




FIG. 4: Sketch defining the distance from the interface s. For 
crystalline particles, the distance s is found to the closest non- 
crystalline particle (any particle with iV < 8). This distance 
is defined to be negative. For liquid particles, the distance 
s is the closest distance to the nearest crystalline particle 
(N > 8). Based on this definition, the region — d < s < d is 
excluded, as no two particles are closer than their diameters d. 
Thus, we shift the distance a by +d/2 for crystalline particles, 
and by — d/2 for liquid particles. This leaves a gap of size d, 
reflecting that the interface truly lies between the crystalline 
and liquid particles. 



sharper, with Win-np = 6. Id; this sharpness agrees with 
previous results [2|, |6( ■ 

However, Fig. [T] implies that measurements that aver- 
age over x and y will artificially broaden the interface, 
given that the interface is not flat. Due to the fractal 
structure of the capillary waves, the interfacial width de- 
pends on the total system size. Our measured widths de- 
pend on the size of the observation region, and would in- 
crease if the observed region was larger and encompassed 
more of the sample [7] . Our data allow us a closer exam- 
ination of the detailed behavior near the interface. For 
each particle, we calculate the distance s to the instanta- 
neous interface position using the procedure illustrated 
in Fig. |H 

To examine the transition from order to disorder, we 
calculate the average number of ordered neighbors (N ) 
now as a function of s, and plot the result in Fig.[5Ja). By 
definition, N > 8 for crystalline particles and N < 7 for 
liquid-like particles, thus producing a gap in (N ) at the 
interface. Figure [5ja) shows that based on the ordering, 
the interface is quite sharp, with a width W10-90 = 1.5d, 
in contrast to the interpretation of Fig. [3] By this mea- 
sure, the interface is locally a sharp, well-defined inter- 
face, as has been conjectured 0, 0|- 

Figure E{b) shows the mobility Ar/d, measured across 
the interface (defined for two different time scales At). 
Here we see the interfacial width is VF10-90 = 2.4d, sig- 
nificantly broader than that determined by structure but 
still fairly narrow, supporting the interpretation of a me- 
andering but sharp interface. On the liquid side, the mo- 
bility is still increasing for s > d where (N ) has reached 
the steady-state liquid value; the halfway point of the 
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FIG. 5: Three views of the intrinsic interface profile, (a) 
Bond-order profile (solid line), showing the average number 
of ordered neighbors of particles as a function of s/d. s 
is the distance from the interface, as described in the text 
and Fig. [4] The dotted line is a hyperbolic tangent fit with 
W^io-90 = 1.5d. (b) Average distance Ar moved in a given 
At as indicated, as a function of s/d. The dotted lines are 
hyperbolic tangent fits, each with Wio-90 = 2 Ad. Here s is 
the distance measured at the particle's initial position at the 
start of the displacement Ar. (c) Nondimensional number 
density d 3 n as a function of s/d. 



hyperbolic tangent fit is at s = 0.2d for the At = 75 s 
data and at s = 0.8d for the At = 750 s data. This result 
is understandable as the crystalline particles form a rela- 
tively immobile region: thus the more mobile liquid-like 
particles are slowed by proximity to the crystal, similar 
to a hydrodynamic "no-slip" condition. Examining the 
individual components of motion (Ax, Ay, Az) we find 
no significant differences, indicating that the motion is 
isotropic, as seen in simulations [2j. While our data for 
At > 750 s are noisy, we can estimate that the diffusiv- 
ity on the liquid side is approximately 50 times smaller 
than the diffusivity in a dilute suspension, in reasonable 
agreement with hard sphere simulations [2|. 

To examine the intrinsic density profile, we plot the 
number density d 3 n(s) in Fig.OJc). The crystal structure 
has a fixed orientation, and thus relative to the curving 
interface, the layers seen in Fig. [3ja) diminish away from 
the interface (the region s < 0). On the liquid side (s > 
0), layering relative to the intrinsic interface is seen with 
two clear layers, fewer layers than the crystal side. This 
reinforces that the layering shown in Fig. 03a) for the 
region z > is mostly due to the uneven interface. 



DISCUSSION 

Note that in this work we have studied slightly charged 
particles rather than ideal hard spheres. Perhaps sur- 
prisingly, our results compare reasonably well with hard 
sphere simulations such as Refs. [2j, la |3, |24j . Thus the in- 
fluence of the charges shift the phase boundaries as noted 
above (a freezing volume fraction ^freeze = 0.49 for our 
experiments as compared to 0.545 for hard spheres), but 
the charges do not seem to strongly modify the structure 
or dynamics of the interface. It is likely that particle mo- 
tion is somewhat slowed due to smaller effective cage sizes 
[20 ] and thus modifying the plateau levels of Fig. \5^b) ■ 
However, this does not seem to modify the extent of the 
transition region of the interface, and indeed one might 
expect that a long-range interaction would only broaden 
the interface, thus making it more striking that we ob- 
serve such narrow interface widths. 

We have studied the equilibrium crystal/liquid inter- 
face of a dense colloidal suspension. We observe a rough 
interface with width 6-10c? depending on the property ex- 
amined. Capillary waves account for some of this broad- 
ening, although overall the mean square fluctuations from 
capillary waves have magnitude gh(At = 0) ~ l.lrf 2 . 
Thus, there may be a base structure comprised of crys- 
talline facets, as suggested by the darkest particles shown 
in Fig.QJb). These facets may be slower to rearrange, and 
their potential existence likely broadens the observed in- 
terface width. 

We also studied the profile in the direction locally per- 
pendicular to the rough interface, and find that this in- 
trinsic interface is much sharper, 1.5c? and 2 Ad based on 
structure and dynamics, respectively. These results em- 
phasize that the observed broadness of the interface is 
only apparent, and depends on the measurement length 
scales [3] and time scales @] . The underlying spatial tran- 
sition from one phase to the other is quite sharp, confirm- 
ing the classic picture [H, 0] ■ Our data provide a useful 
test for future models of the intrinsic profile. 



MATERIALS AND METHODS 

Our samples are composed of colloidal poly- 
methyl(methacrylate)_particles, sterically stabilized to 
prevent aggregation [13j |. The particles have mean di- 
ameter d — 2.30 ± 0.02 /jm and a polydispersity of ap- 
proximately 5%. The solvent is a mixture of cyclohexyl- 
bromide and decalin, chosen to closely match the density 
and index of refraction of the particles, with viscosity 
770 = 2.25 mPa-s 2JJ. However, we intentionally slightly 
density mis-match the fluid by adding an excess of decalin 
(which is more dense than the particles). This lets us 
gravitationally induce a crystal-liquid transition within 
our sample chamber. Our particle/fluid density mis- 
match is less than Ap s=s 5 x 10 g/ml, giving us a grav- 
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itational scale height ksT/ ApVg > 130/im, using the 
temperature T (295 K for these experiments), the parti- 
cle volume V = f-d 3 , and the gravitational acceleration g. 
In our solvent, the particles have a slight charge, shifting 
the freezing transition volume fractions to <^f rC czc = 0.43 
and the melting transition to mc it = 0.49, similar to 
values for other weakly charged colloidal PMMA experi- 
ments [HI . The existence of the charge is also apparent 
from comparisons of the pair correlation function g(r) 
taken at different volume fractions. For hard spheres, 
the first peak of g(r) should always be at r = d, whereas 
for our samples we find the peak at r > d and the peak 
position varies slightly with <j). However, this peak posi- 
tion is still fairly close to d; for example, in the crystalline 
region we find the peak at r = 1.04d and in the liquid 
region we find the peak at r = 1.15d, suggesting that the 
charging is not too extreme 28, 29[. 

Our microscope sample chambers have dimensions 0.4 
mm x 3 mm x 30 mm, and are stored with the long di- 
mension oriented vertically. Because of the slight density 
mismatch between the solvent and the particles, grav- 
ity sets up a slight concentration gradient. We let the 
samples equilibrate in the sample chamber for more than 
one month before taking data. We make one important 
modification to our microscope: the microscope base is 
tilted 90° so that the objective lens points horizontally 
[3(ij |. The sample is placed on the microscope stage with 
the orientation kept the same as the storage conditions. 
Thus we study a completely equilibrated sample at a sta- 
ble crystal-liquid interface, and the direction of gravity 
(— z) points perpendicular to the optical axis. 

The particles are dyed with rhodamine 6G, so that 
they can be viewed with a laser scanning confocal micro- 
scope ljl 27 1 . We acquire images of size 20 x 55 x 60 /xm 3 
once every 25 s, where the long direction (z) is parallel to 
gravity. The images were taken from 20 to 40 /jm away 
from the coverslip, to avoid direct influence of the walls. 
As noted above, the full sample size is 0.4 x 3 x 30 mm 3 , 
where the smallest dimension is x and the largest direc- 
tion is z; thus the imaged region is only a small volume 
within the sample. The 3D images are analyzed to de- 
termine particle positions with a resolution of 0.1 /im 
parallel to the optical axis (x) [2?], [3l[ . and a resolution 
of 0.05 /jm perpendicular to the optical axis (y and z) 
Because this is a dense sample, particles do not move far 
between images, and thus we follow their motion using 
standard particle tracking techniques [3lj ] . 
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